Lagrangian Investigation of Two-Dimensional Decaying Turbulence 
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We present a numerical investigation of two-dimensional decaying turbulence in the Lagrangian 
framework. Focusing on single particle statistics, we investigate Lagrangian trajectories in a freely 
evolving turbulent velocity field. The dynamical evolution of the tracer particles is strongly domi- 
nated by the emergence and evolution of coherent structures. For a statistical analysis we focus on 
the Lagrangian acceleration as a central quantity. For more geometrical aspects we investigate the 
curvature along the trajectories. We find strong signatures for self-similar universal behavior. 
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Introduction Beyond its importance for geophysi- 
cal and astrophysical applications two-dimensional tur- 
bulence often serves as a paradigm of a complex, self- 
organizing system. Despite its spatio-temporal complex- 
ity freely decaying two-dimensional turbulence exhibits 
a stunning degree of coherence. The dynamical evolu- 
tion of the turbulent vorticity field is dominated by the 
emergence, coalescence and nonlinear interaction of cir- 
cular and spiral-like vortices (see e.g. [l] and references 
therein) . 

During the decay process these coherent structures 
tend to organize in a self-similar way 0, H, [l] j for ex- 
ample the temporal evolution of their density, the mean 
absolute value of the circulation and the mean vortex 
radius follow power-laws [J, indicating that the decay 
process might display some degree of universality. 

The Lagrangian frame of reference offers a natural ac- 
cess to the description of turbulent flows. However, re- 
sults on the Lagrangian description of two-dimensional 
turbulence are sparse. Hence we are led to study the 
evolution of tracer particles in the two-dimensional de- 
caying system. We find, as expected, that their evolution 
is strongly dominated by the coherent vortices. 

Regarding a statistical analysis we focus on the accel- 
eration as well as on the curvature along the trajectories. 
We find strong evidence for a self-similar temporal evo- 
lution of some of the quantities. 

The system under consideration is governed by the 
nonlinear effects of vortex interaction and merging as well 
as by viscous forces. Thus please note that our results 
cannot be referred to in the framework of the universal 
decay theory proposed by Carnevale et al. 0, Q as we 
operate in a different parameter range, where viscosity 
cannot be neglected. 

The remainder of this article is structured as follows. 
After summarizing details on numerical issues and simu- 
lation parameters we procede to a qualitative discussion 
of trajectories. We then present statistical results of the 
Lagrangian acceleration and curvature. 
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Figure 1: Temporal evolution of the energy spectrum for 
t £ [0,16]. High modes get damped out while the energy 
contained in the low modes increases. 




'Electronic address: ' mwilczek@uni-muenster.del 



Figure 2; Vorticity field for t — 0.3 and t — 4. The number of 
vortices decreases due to vortex coalescence and destruction. 



Numerics and Simulation Details We solve the 
two-dimensional vorticity equation 

^{x, t) + {u{x, t) ■ V)uj{x, t) = vl\uo{x, t) (1) 
at 

by means of a standard pseudospectral code on a 1024^ 
grid with periodic boundary conditions. Table [J sums 
up the major simulation parameters. We suppress alias- 
ing errors by a spherical mode truncation according to 
Orszag's famous two thirds rule. To obtain reliable sta- 
tistical results we immerse approximately 16000 tracer 
particles into the flow and follow them throughout the 
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L U Ret dt T Mra iVrc 

1024^ 2-K 0.0004 12500 0.0005 25 16384 40 

Table I: Major simulation parameters. Number of grid points 
A^'^, box length L, viscosity v, initial Reynolds number Rei — 
^^L-iEsis ^ timestep dt, total simulation time T, number of tracer 
particles per realization A'tra and total number of realizations 

iVre. 



decay process. The Eulerian fields are interpolated by 
a bicubic scheme. The acceleration is explicitly calcu- 
lated each timestep by evaluating the right-hand side of 
the Navier-Stokes equation. The resulting field is then 
interpolated at the particle positions. 

In order to increase statistical quality we additionally 
perform an ensemble average over forty independent re- 
alizations of the flow (about 19 gigabytes of Lagrangian 
tracer data, 0(1600) hours of cpu time). The initial con- 
ditions are obtained from a forced turbulence simulation 
that has reached statistical equilibrium in the following 
manner. A master field from such a forced simulation 
is taken and forty copies are made. For each copy a 
forced simulation with a duration of about five large eddy 
turnover times is performed with a differently initialized 
random number generator for the forcing (see [6] for de- 
tails on the forcing). We thereby gain forty statistically 
independent initial conditions for the subsequent decay- 
ing turbulence runs. We would like to stress that by the 
large amount of ensemble realizations we achieve an ex- 
traordinarily high statistical quality. One additional run 
with doubled total simulation time was performed for 
checking and visualization purposes. For a more detailed 
description of the numerical issues see reference @. 

The initial conditions exhibit a clear inertial range ac- 
cording to Kolmogorov's predictions as depicted in figure 
[T] (here, no ensemble average was performed). The time- 
resolved energy spectrum reveals that the initial inertial 
range scaling is quickly destroyed as the high modes get 
damped out. Simultaneously energy is transported into 
the low modes by nonlinear mode interaction of the in- 
verse energy cascade. In physical space this is mirrored 
by the decay of small-scaled structures and the emergence 
of bigger coherent vortices as depicted in figure [2] 

Trajectories Turning now to a Lagrangian descrip- 
tion of the turbulent field, figure [3] shows some sample 
trajectories. As characteristic for a turbulent chaotic sys- 
tem two initially neighboring trajectories ((1) and (2)) 
separate quickly as they get caught into the velocity field 
of different vortices. The spiraling motion of the parti- 
cles indicates that vortex trapping is generic for particles 
in two-dimensional decaying turbulence. As the vortices 
merge into bigger structures the radius of the spiral mo- 
tion increases on average over time. We postpone a more 
quantitative discussion of this point to the following sec- 
tion. 

While the chaotic advection can significantly increase 
the pair separation also the opposite takes places. For 
example trajectory (1) and trajectory (3) subsequently 
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Figure 3: Typical particle paths in decaying turbulence. The 
spiraling motion indicates vortex trapping events. The num- 
bers denote the initial position of the three particles. 
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Figure 4: Time-resolved trajectory in decaying turbulence. 
Vortex trapping results in an oscillating behavior of the com- 
ponents of position, velocity and acceleration. 



enter the velocity fields of the same vortex leading to 
a sticking effect; the vortex trapping event causes these 
two particles to stay close for a comparably long time 
although they started from far separated places in the 
turbulent field. 

Figure |4] shows the time-resolved trajectory (1). Af- 
ter the initial emergence of medium-scaled vortices from 
the rather rough and unstructured initial vorticity all the 
components of space, velocity and acceleration show the 
signatures of vortex trapping by their oscillating behav- 
ior. The strength and also the frequency of the oscilla- 
tions decrease throughout the decay process as the co- 
herent structures increase in size and the velocity field 
gets smoother and damped because of viscous effects. 
All these observations show that Lagrangian particle dy- 
namics and hence statistics is strongly dominated by the 
motion in single coherent structures. 
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Figure 5: Time-resolved acceleration pdf. The upper fig- 
ure shows the acceleration pdf's for different times vertically 
shifted. The lower figure shows the same pdf's normalized by 
their standard deviation resulting in a collapse of the curves. 
Time-windows over which averaging is performed are indi- 
cated. 



Statistical Results After this rather qualitative dis- 
cussion let us now turn to the central statistical observ- 
able of this work. The Lagrangian acceleration of a par- 
ticle initially located at y is defined by 

a{t, y) = X{t, y) = [-Vp{x, t) + v^u{x, t)]^^x(t y) ■ 

For a more detailed dynamical description this accel- 
eration can be decomposed into a acceleration paral- 
lel to the current velocity and a acceleration perpen- 
dicular to the current velocity. With ity = u/\u\ and 
u±_ = {—Uy,Ux)/\u\ the decomposition reads 



a{t,y) 
thereby defining 



0'\\it,y)u,i + aj_{t,y)uj_, 



a|| = a ■ M|| 



and 



= a ■ uj_. 



(3) 



(4) 



Let us start with the time-resolved pdf /(a^y) of a single 
component of the acceleration, which is depicted in figure 
[51 Owing to the statistical isotropy of the flow we improve 
our statistics by averaging over both spatial directions. 



Figure 6: Time-resolved moments of the acceleration pdf. The 
moments clearly show power-law behavior. The inset shows 
a comparison with power-law functions. The algebraic expo- 
nents are displayed. 



The time-resolved acceleration pdf's clearly exhibit large 
deviations from a Gaussian distribution as expected from 
the intermittent nature of turbulent acceleration. They 
exhibit almost exponential tails. However, when normal- 
ized by their standard deviation all of the curves col- 
lapse over the selected time interval thereby indicating 
a universal temporal evolution of the pdf. This view is 
supported by figure [5] which shows the time-resolved mo- 
ments of the acceleration pdf. By inspection of this figure 
the scaling regimes roughly can be estimated as the time 
interval for t S [0.1,10]. Note that this choice is rather 
conservative. The size of this time interval will surely 
depend on physical parameters like viscosity or initial 
Reynolds number. Additionally statistical issues might 
matter; as the turbulence decays the number of vortices 
rapidly drops. Consequently there are fewer and fewer 
coherent structures contributing to the overall statistics 
resulting in a decreasing statistical quality. 

Figure [6] also shows a comparison with power-law func- 
tions {a"y) t''" with the exponents C„. The numerical 
values roughly suggest (^2n nC,2 for the selected mo- 
ments as necessarily required for a self-similar evolution 
of the pdf's. Please note that we do not claim these ex- 
ponents to be universal for all turbulent flows as they 
should dependend on the Reynolds number, the type of 
numerical viscosity applied in the simulation and on the 
exact initial condition. For a detailed investigation of the 
impact of different initial conditions on the decay process 
see 0. 

We now proceed to a discussion of the parallel and 
perpendicular component of the acceleration. Figure [7] 
shows a comparison of the pdf of the acceleration com- 
ponents a 



xy, 



and a±_ . One can see that the perpendic- 
ular component is slightly more intermittent than the xy- 
averaged component. However, the parallel component 
is clearly less intermittent. This difference remains when 
the pdf's are rescaled to their standard deviation (not 
depicted) indicating that the three components have a 
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Figure 7: Pdf's of the acceleration components a^^y, a\\ and 
a± for the time-window 4-6. In comparison with a^y the per- 
pendicular component is more intermittent while the parallel 
component is less intermittent. 
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Figure 8: Temporal scaling for the second and fourth moment 
of the acceleration components a^y, a|| and a±. While the 
scaling exponents seem to be almost identical the moments 
differ by a multiplicative factor. 



fundamentally differing functional form. The physical in- 
terpretation of these results is quite straightforward. The 
fact that the perpendicular component of the acceleration 
is more intermittent than the parallel one shows that the 
spiraling motion in the vortices is a major contributor 
to the intermittent nature of the turbulent acceleration. 
Hence this is statistical evidence for the dynamical im- 
portance of vortices regarding Lagrangian dynamics. 

We checked that also the pdf's of the components a\\ 
and aj_ collapse over the selected interval when normal- 
ized by their standard deviation (not depicted) and hence 
also the moments show a power-law behavior. This is 
exemplified in figure [8l where the second and fourth mo- 
ment of the corresponding pdf's are shown. This inspec- 
tion reveals that the moments show nearly identical tem- 
poral scaling behavior, i.e. (a"^) = Cxyf^", (ap = c\\ i^" 
and (a") = c±t'^" with identical C„, but differ by their 



Figure 9: Log-log plot of the mean curvature as a function of 
time. The increasing spiral radii of the trajectories result in 
a decaying curvature. 
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Figure 10: Log-log-plot of the time-resolved curvature pdf. 
The pdf's show an algebraic decay with an exponent close to 
-2.25. 



prefactors Cxy, cy and c±. That means the differing func- 
tional form can be traced back to a differing weighting 
for each moment. The emerging picture is quite inter- 
esting: while the functional form of each of these pdf's 
seems to be determined by an initial weighting factor for 
each moment, their temporal behavior is universal. 

To close this section let us turn to more geometrical 
properties of the trajectories. We investigate the curva- 
ture of the trajectories defined by 



(5) 



Figure [5] shows the mean curvature, averaged over all 
trajectories. Although this quantity is extremely fluctu- 
ating, that is because of the l/v? dependence, the aver- 
age trend of a decreasing curvature is clearly visible. As 
the curvature is proportional to the inverse of the spi- 
ral radius, the emergence of bigger and bigger coherent 
structures therefore on average leads to increasing spiral 
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radii of the trajectories. This is the statistical approval 
for the qualitative discussion of the trajectories in the 
preceding section. The time-resolved pdf of the curva- 
ture contains even more information and is shown in fig- 
ure [TOl The time- resolved pdf's arc fairly similar to each 
other. However, the probability for a high curvature (a 
small spiral radius) decreases over time, as expected for 
coalescing and growing vortices. Interestingly, all of the 
pdf's show an extraordinarily clear algebraic decay with 
an exponent close to —2.25. 

Summary To sum up, we presented a detailed inves- 
tigation of two-dimensional decaying turbulence in the 
Lagrangian frame of reference. Our results reveal that 
the Lagrangian dynamics in this system is strongly influ- 
enced by long periods of trapping in coherent structures. 
The spiraling motion of the particles is a consequence. 

We focused on the Lagrangian acceleration as a central 



statistical quantity finding that the turbulent fields tend 
to organize in a self-similar way. Hence we are able to 
identify a time interval where a self-similar scaling of the 
acceleration pdf's and corresponding moments holds. A 
decomposition of the acceleration into components paral- 
lel and perpendicular to the velocity reveals that the in- 
termittent nature of the acceleration is partly due to the 
centripetal accelerations caused by the coherent struc- 
tures. Closely connected to the acceleration is the curva- 
ture of the trajectory, adding a more geometrical point 
of view. We found the mean curvature decaying alge- 
braically in time, consistent with the qualitative picture 
drawn in the first sections. The time-resolved pdf's of 
the curvature show a clear power-law behavior. 

Our investigation gives evidence that two-dimensional 
decaying Lagrangian turbulence is dominated by self- 
similar scaling laws, even when viscosity is not negligible. 
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